Tunneling measurement of quantum spin oscillations 
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We consider the problem of tunneling between two leads via a localized spin 1/2 or any other microscopic 
system (e.g., a quantum dot) which can be modeled by a two-level Hamiltonian. We assume that a constant 
magnetic field Bo acts on the spin, that electrons in the leads are in a voltage driven thermal equilibrium and 
that the tunneling electrons are coupled to the spin through exchange and spin-orbit interactions. Using the 
non-equilibrium Keldysh formalism we find the dependence of the spin-spin and current-current correlation 
functions on the applied voltage between leads V, temperature T, Bo, and on the degree and orientation m Q 
of spin polarization of the electrons in the right (a =R) and left (a =L) leads. We show that a) The spin-spin 
correlation function exhibits a peak at the Larmor frequency, ojl, corresponding to the effective magnetic field 
B acting upon the spin as determined by Bo and the exchange field induced by tunneling of spin polarized 
electrons, b) If the m a 's are not parallel to B the second order derivative of the average tunneling current 
I(V) with respect to V is proportional to the spectral density of the spin-spin correlation function, i.e., exhibits 
a peak at the voltage V — hujh/e, c) In the same situation when V > B the current-current correlation 
function exhibits a peak at the same frequency, d) The signal-to-noise (shot noise) ratio R for this peak reaches 
a maximum value of order unity, R < 4, at large V when the spin is decoupled from the environment and the 
electrons in both leads are fully polarized in the direction perpendicular to B, and e) R <C 1 if the electrons are 
weakly polarized, or if they are polarized in a direction close to Bo, or if the spin interacts with the environment 
stronger than with the tunneling electrons. Our results of a full quantum-mechanical treatment of the tunneling- 
via-spin model when F > B are in agreement with those previously obtained in the quasi-classical approach. 
We discuss also the experimental results observed using STM dynamic probes of the localized spin. 
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PACS numbers: 03.65.Xp, 03.65.Ta, 73.40.Gk, 73.63.Kv 



I. INTRODUCTION 

Interest in quantum information processing has brought sig- 
nificant attention to the problem of measurement of tunneling 
currents via a microscopic system that can be modeled by a 
two-level Hamiltonian (such as a quantum dot, or a molecule 
or an atom with a localized spin HI MlJ). In the case of 
a single spin, measurements of the tunneling current in such 
a system provide information on spin orientation and its dy- 
namics, and they constitute an example of indirect-continuous 
quantum measurement jjj] . A fundamental question that arises 
is what signatures of the spin dynamics are encoded in the tun- 
neling current and how this current affects the spin dynamics. 

Scanning tunneling microscopy (STM) experiments j^l 0] 
on a single molecule with a spin, in the presence of a mag- 
netic field Bo, have reported a peak in the current noise power 
spectrum (i.e., current-current correlation function) V(u>) at 
the Larmor frequency lol = jBq, where 7 is the gyromag- 
netic ratio. Experiments were done at room temperature and 
the authors found that the signal-to-noise ratio R (ratio of the 
power at the peak frequency to the shot noise power) exceeded 
unity and was almost independent of the orientation of the 
applied magnetic field Bo- In the non-relativistic approach 
the tunneling electrons couple to the spin by exchange inter- 
action. In this case, electrons with a spin polarization along 
Bo do not couple with the oscillatory components of the spin 
(which are perpendicular to Bo). The experimental results 
0, 0] are difficult to explain in the framework of a single- 
spin non-relativistic model since electrons in the leads were 
polarized by the same magnetic field which acted on the spin. 
Possible relevance of spin-orbit interaction to explain these 



data was discussed by Shachal and Manassen and later on b 
Balatsky and Martin 01 . Recently, Levitov and Rashba 
noticed that in systems with low space symmetry (as in dots 
or a molecule near a surface) the nonvanishing orbital moment 
of electrons near the spin (in the dot or molecule) provides a 
strong coupling of the tunneling electrons to the spin, and they 
speculated that this mechanism may lead to a significant effect 
of the spin oscillatory component on the tunneling current. In 
fact, in order to fit the experimental data |5l |6[ one needs a 
model of electron tunneling via a single spin which is able to 
explain not only the existence of a peak in the current power 
spectrum, but also a significant signal-to-noise ratio R > 1 
almost independent of the orientation of Bo- 

In the following we consider a model with exchange (non- 
relativistic) coupling of a single spin 1/2 and the tunneling 
electrons, including spin-orbit coupling. In the framework of 
this model we analyze the dependence of R and linewidth T 
on the applied voltage V between leads, the applied magnetic 
field Bo, the temperature T, and on the degree and orienta- 
tion m a o of electron polarization in the right (a =R) and left 
(a =L) leads in the steady state. This state establishes in a 
transient time after the voltage or tunneling matrix elements 
are switched on. 

The problem of tunneling via a single localized spin 1/2 
is similar to the problem of tunneling via a double quantum 
dot (two-level system) extensively studied by Korotkov and 
Averin |3| and by Ruskov and Korotkov U Indeed, they 
were interested in the question of how quantum oscillations 
in the two-level system, representing a qubit, may be detected 
by a tunneling current. Such a system can be described by 
a model of electron tunneling via a spin when leads elec- 
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trons are fully polarized. Using the Bloch equations describ- 
ing the ensemble-averaged evolution of the density matrix for 
the coupled qubit-detector system, they obtained oscillations 
at the Larmor frequency in the noise power spectrum V(u>), 
with R < 4 in the case of weak interaction between the qubit 
and the detector. A similar quasi-classical approach for tun- 
neling of electrons via a single spin was used in Ref. Q|. In 
such an approximation, the action of the electrons on the spin 
is replaced by the action of an effective classical magnetic 
field with a shot noise spectrum depending on the tunneling 
current between the leads. This approach has sense at high 
currents (voltages), while at low currents it does not account 
for the quantum nature of the tunneling electrons and hence 
a more elaborate solution of the quantum transport equation 
for both spin and electrons is desired. Here the situation is 
formally similar to the case of the photo-electric effect where 
in some circumstances a full quantum mechanical treatment of 
the problem provides the same information as a semi-classical 
approach, while in other cases, such as the measurement of the 
statistics of the photo-current fluctuations, the predictions of 
the two approaches differ qualitatively . 

An important step toward a full quantum mechanical treat- 
ment of electron tunneling via a spin was recently given by 
Parcollet and Hooley |2]. Assuming that the spin interacts 
only with the tunneling electrons, and electrons in the leads 
are in a voltage driven thermal equilibrium, they computed the 
spin magnetization of a quantum dot in the two-leads Kondo 
model as a function of the temperature of the electrons, mag- 
netic field, and voltage between the leads (which exceeds the 
Kondo temperature Tk) in the steady state, i.e., in the long 
time limit. They considered the situation where electrons in 
the leads are polarized in the same direction as the magnetic 
field acting on the spin, rriRo = hilo = mo, mo 1 1 Bq- Their 
calculations confirmed previously obtained quasi-classical re- 
sult fioll that, in the steady state even at zero order in the spin- 
leads coupling, the magnetization is not given by the thermal 
equilibrium expression A/ eq = (1/2) tanh(_Bo/2T) when the 
spin is more strongly coupled to the leads than to any other 
thermal bath. At long times, the state of the spin is completely 
determined by the characteristics of the tunneling electrons, 
and since the system is out of equilibrium, the steady state 
thus achieved is described by a steady state distribution func- 
tion which differs from the Gibbs distribution. It turns out 
that this distribution function depends upon the voltage be- 
tween the leads. They also showed that the correct way to 
calculate the perturbative corrections to the spin distribution 
function and the spin decoherence rate T in the steady state 
is by solving a quantum transport equation self-consistently. 
We will extend the treatment of Parcollet and Hooley to arbi- 
trary orientations of mjo, and hilo with respect to Bo, and 
account for a direct tunneling of electrons between leads. 

The main goals in this article are (i) to examine un- 
der what tunneling conditions the current-current correlation 
function V(lS) exhibits a peak at the Larmor frequency in 
the steady state, and (ii) to calculate r(V, T, Bo, m Q o) and 
R(V, T, Bo, m Q o) in the regime where eV ^> Tk, with e de- 
noting the charge of the electron. As in Ref. |01 we use the 
Keldysh formalism JTllfl2il and the Majorana-fermion repre- 



sentation fl3ll for the spin to find the spin distribution function 
and the current-current correlation function. 

We note, from a more general perspective, that the prob- 
lem of electron tunneling via a spin represents an example of 
indirect quantum measurement. The spin is probed by the tun- 
neling electrons whose correlation function is measured in a 
continuous fashion with a classical apparatus (which is not in- 
cluded). The goal of quantum measurement is to determine 
the frequency of precession of the isolated spin (energy sep- 
aration in the isolated two-level system) and to obtain infor- 
mation on the initial state of the system. However, when the 
spin is decoupled from the environment, tunneling via the spin 
changes its state as the measurement goes on. After switch- 
ing on the voltage, or the tunneling matrix elements, there is 
an initial period of time where the tunneling current depends 
upon the initial state of the spin. After some transient time 
the steady state, which is in general independent of the ini- 
tial conditions, establishes. In the steady state the frequency 
of precession is renormalized by the tunneling electrons, and 
the width of the peak (at the precession frequency) and cor- 
responding signal-to-noise ratio depend upon the precession 
frequency and characteristics of the tunneling electrons. To 
extract information about the precession of the isolated spin 
one needs to have complete information on the tunneling elec- 
trons. We show how to relate the precession frequency of the 
isolated spin to results of current measurements in the steady 
state (peak frequency, width of the peak and signal-to-noise 
ratio). To obtain information on the initial state of the spin, 
measurements of the current during the initial transient period 
are needed. Here we do not consider the spin dynamics in this 
short time interval. 

The plan of the article is the following. In the next section 
we present the model for tunneling via a single spin 1/2, with 
and without spin-orbit coupling, and a model describing tun- 
neling via a two-level quantum dot. Then we introduce Majo- 
rana fermions and the Keldysh technique. In the subsequent 
section we will compute the spin distribution function, broad- 
ening of the spin precession due to the tunneling of electrons, 
and the spin-spin correlation function. Then we will calculate 
the dependence of the tunneling current on the spin and deter- 
mine the current-current correlation function V{uj). Finally, 
we discuss the STM experimental results and compare our re- 
sults to those obtained in the quasiclassical approach 01 00] 
and other theoretical works QBCal- 



II. THE MODEL 

For the system consisting of a spin coupled to the leads 
by an exchange mechanism (in the non-relativistic approxi- 
mation, i.e., neglecting the spin-orbit interaction), we use the 
Hamiltonian of the two-leads Kondo model (jl 01 with a di- 
rect tunneling term included (which we call tunneling-via- 
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spin (TvS) model) 

Tt = 7i e + TL S + H.T , H-T = 7~iref + T~ttn (1) 

W e = y ] [ e na^crcr' — ' IT ™ , l C Lii c aniT'! 

W s = — g/isBo • S, 

Wref = ^ ] c <Lncr {T aa ) aa r C an i a i , T aa — T^ Q ^ S ■ CTaa' j 
(Trl)ctct' = T^aa' + Tjy* S ■ <T<ra< , 

where c^^. (c Qn(T ) creates (annihilates) an electron in the left 
or right lead (depending on a £ {L, R}) in the eigenstate 
71, and with spin a. Further, e na = e n — p Q , where e„ is 
the energy in the state n and p a is the chemical potential 
in the lead a, while er represents the three Pauli matrices. 

Ill • -^rr an ^ ^lr are tunneling matrix elements due to 
the exchange interaction for the electron tunneling from the 
leads to the molecule with the spin 1/2, while To is the di- 
rect tunneling matrix element. We take them as real numbers. 
The spin localized in the molecule is described by the oper- 
ator S = (S x ,S y ,S z ). Figure Q sketches the physical setup 
we want to study and which basically represents the model 
Hamiltonian H. 

In the following we use energy units for the bias voltage V, 
i.e., we denote eV by V. Also, by B we mean gpsT, write 
T instead of fc^T, and u represents hw. 

We describe the leads by a free electron gas with a density 
of states p(e) of bandwidth D, so that e n is the bare energy of 
the electron in the eigenstate n, the same in both leads. We as- 
sume weak tunneling, Tq pg, \ 2 p\\ <C 1, where po is the 
density of states per spin (DOS) of the leads at the Fermi level 
(when the leads are different p\ = p$p^ where pfj is me DOS 
in the lead a). The applied voltage is represented by the differ- 
ence in chemical potentials, pl — PR = V. We also assume 
V » T K « \T^\eM-\T^ ] r 2 Po 2 \ and T » T$ 
(a, (3 £ {R, L}). We assume in the following that Bpo <C 1 
and Vpo <C 1; in other words, the lead-electron bandwidth 
D ~ 1/po is the largest energy scale. The leads are supposed 
to be in a voltage driven thermal equilibrium at temperature 
T. Href describes spin-flip scattering of an electron from one 
lead back into the same lead, and H tr represents both direct 
and spin-assisted tunneling between leads. In the case of fully 
polarized leads, the summation in the Hamiltonian is taken 
only over one type of electron spin. 

In the model with spin-orbit interaction proposed by Levi- 
tov and Rashba |8] the spin precession induces electron den- 
sity oscillations inside the dot, and thus an electric field inside 
and around the dot, E = (n A (S A L)), where L is the orbital 
momentum, while n is a polar vector allowed by the symme- 
try of the system. This electric field modulates the tunneling 
barrier between the leads and the dot. Hence, the tunneling 
matrix element acquires, in addition to the non-relativistic ex- 
change term, a dependence on the localized spin S which is 
similar for both spin components a of the tunneling electrons. 




FIG. 1: Schematics of the physical systems represented by the 
Hamiltonian Ti, and measurement process involved in determining 
the current (I(t)) and current-current correlation function (I(t)/(0)) 
which contain encoded information about the spin S. The electronic 
tunneling current is established by a dc voltage V. 



To describe such relativistic corrections we need to add to the 
transfer matrix elements in the Hamiltonian, Eq. Q, the spin- 
orbit terms. The compound matrix elements are finally given 

as 

(Tqlol)(7(7' S ■ (T aa > -\- T^ a ^S • 1 b aa ' , (2) 

(T RL ) CTCT , = (T + T^ o) S • 1) <W + T£ L x) S ■ <w. (3) 

Here the unit pseudo-vector 1 indicates what projections of the 
localized spin are modulating the tunneling matrix elements 
depending upon the geometry of the system. We can estimate 
~ r 2 Eo, where To is an energy of the order of an atomic 
energy, while r is a small dimensionless parameter which de- 
scribes the weak overlap of the lead-electron wave-function 
to that of the molecule with the spin. For spin-orbit tunneling 
we have an additional small relativistic multiplicative factor 
(3 2 = v 2 /c 2 , i.e., T£p ~ /3 2 T^, where v is a typical elec- 
tron velocity. 

Finally, we presume that the spin also couples to the envi- 
ronment due to the spin-phonon coupling, the coupling with 
nuclear spins, etc. To account for this spin relaxation mecha- 
nism we introduce the relaxation rate r env . 

In order to put our model in perspective let us compare it 
to other models used in the literature. For the quantum dot 
studied in Refs. |3l |4[ the Hamiltonian for spinless fermion 
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tunneling is given by 



(4) 



b 



-(4d 2 +4^l) + ff tr, 



n tr = J2i M + -^(4* - 4*)] 4™%™ + H - c -> 

where d\ 2 (g^ 2 ) creates (annihilates) an electron in the state 
1 or 2 of the dot. Here b z /2 is the energy splitting between 
states 1 and 2, b x /2 is the tunneling matrix element connecting 
them, and tunneling of electrons between the leads depends on 
the population of the states 1 and 2. The difference between 
this and the TvS model (with H re f=0) is that the couplings to 
the leads only depend upon the difference of fillings of states 
1 and 2 or, making correspondence with the spin operators, 
the Hamiltonian H s is 



H s = b x S x + b z S z , 

s x = (i/2)(4d 2 + 4d 1 ), 
s z = (i/2)(4di - 4d 2 ), 



(5) 



while the coupling, Ji tr depends only on one spin component, 
S z of the probed system, instead of all of them as in the ex- 
change term S • cr. Moreover, the leads electrons are fully 
polarized along the z-axis. 



III. MAJORANA-KELDYSH FORMULATION OF 
TUNNELING VIA A SINGLE SPIN 

It is assumed that for sufficiently long times the compos- 
ite (spin-electrons) system reaches a dc (non-equilibrium) 
steady state (which does not depend upon the initial condi- 
tions). We use the Majorana fermion representation for the 
spin and the Keldysh diagrammatic technique lfl2ll to describe 
it. This technique is nowadays a standard method used in non- 
equilibrium phenomena and its basic description can be found 
in standard textbooks such as Ref. fl5ll . Its application to 
Kondo-like problems is described in Ref. Q]. 

Since spin operators are not amenable to the application of 
the Wick's theorem, and hence are not appropriate for dia- 
grammatic techniques, we represent spin- 1/2 operators in 
terms of three Majorana fermions r/^, where [i € {x, y, z}, 
satisfying the following mapping and relations 

S x = -iVyVz, S y = -ir] z r] x , S z = -irj x r) y , (6) 
fe) t = {n^Vv} = <W (7) 

The Hilbert space for Majorana fermions is eight-dimensional 
(8<f) and the transformation from the spin 1/2 to the Majo- 
rana fermion representation may be performed in the follow- 
ing way. In general, the spin expectation values that we need 
to calculate are of the form of a trace, Tr^lS^}], over the 
two-dimensional (2d) spin Hilbert space. Here F is some 
functional of the spin operators. We replace the spin operators 



by Majorana fermion operators using Eq. (|6j and we represent 
the Majorana fermion operators as 



(4 + r)/V2, 



(8) 



where 4 and c„ are the creation and annihilation operators 
of a type fi e {x, y, z} fermion satisfying the canonical anti- 
commutation relation {c^, cj,} = 8^ u . This representation 
guarantees that Eq. Q for the Majorana fermion operators 
are satisfied. A possible orthonormal basis for the 8d Hilbert 
space Hm is <fi = where $0 is the vacuum, ip 2 = 4^o> 
ip 3 = 4$ , ip 4 = 4$ , tp 5 = 44 $0 ' ^ 6 = 4 c z* ' 
(fi7 = 44^o> and ip 8 — c x c y c\^ . % may be writ- 
ten as a direct sum of four orthogonal 2d subspaces, Hm = 
Tti+ © Hi- © H2+ © H2-, which are invariant under the 
action of the spin-1/2 operators S^. Here, basis sets for Hi+, 
Hi-, H 2+ , andH 2 - are: {(ip 6 +icp 7 )/y/2, (ipi + iip 5 ) / V2}, 
{(ipi - icp 5 )/V2, (<p 6 - icp T )/V2}, {(f 2 + 1^3) /V2, (<p 4 + 
i(fg)/V2}, and {(^4 - up&)/y/2, (<p 2 - iip 3 )/V2}, respec- 
tively. In other words, the Majorana fermion representation 
of the spin-1/2 operator is an irreducible one in any of the 
four mentioned subspaces. Furthermore, these representations 
are equivalent to each other. Therefore, we can write down 



1, 



Tr[F{S„}] = -Tr[F{ Vfl }], 



(9) 



where the second trace is computed in the extended Hilbert 
space Hm- (The trace is invariant under a change of basis in 
Hm)- Now one can prove Wick's theorem in the same way as 
it was proved for the free Fermi and Bose fields. 

The time-ordered averages in the Keldysh technique are 
taken with the help of the evolution operator 5c along a closed 
time contour (its part running from — 00 to +00 denoted by + 
and the part from +00 to —00 by — ). One defines the four 
(non-independent) real-time Green's functions 



<V (*,*') = 



-i{m^{t')), 
i(i>\t')i>{t)), 



(10) 



which can be compactly written as a matrix G^(t,t') = 
-i((^{t)^(t'))). Here T and T are the time- and anti-time- 
ordering operators on the + and — parts of the Keldysh con- 
tour and tp represents an arbitrary fermionic field. They are 



related through: G 



g: 



GL , and the brack- 



ets ) have the meaning of either a pure state or a distribution 
of the available phase space of the interacting system. 

The four Green's functions (the "± basis") can be expressed 
in terms of the retarded, advanced and Keldysh Green's func- 
tions ("Larkin-Ovchinnikov (LO) basis") 



G R 



^ (*,*') 



(id 
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by use of the transformation 



r<R r<K 

^■4, 







1 1 

1 -1 



C++ 



G 



4> 



1 1 

-1 1 



In general, the retarded and advanced functions are Hermitian 
conjugates of each other, i.e., G$(t,t') = (G$(t\t))*. Fur- 
thermore, in thermal equilibrium the fluctuation-dissipation 
theorem connects them through 



G$(u)=h eq (u)[G$(w)-G$(uj)], 



(12) 



„ eq (w) 



with the thermal distribution function of fermions /i e _, , 
— tanh(o;/2T). As in Ref. 01 we assume that for sufficiently 
long times the spin reaches a dc (non-equilibrium) steady state 
which does not depend on the initial conditions. In such dc 
steady state the correlation functions depend upon the time 
difference and the new distribution function h(u>) differs from 
h eq (u). It should be determined as a stationary solution of the 
quantum kinetic equation. 

Firstly, we introduce the bare Green functions for electrons 
in the presence of an applied magnetic field (or internal ex- 
change field in the case of magnetically ordered leads) B Q 
which induces a spin polarization of electrons in the direc- 
tions m Q o = B^/Bct- Since the leads are assumed to be 
coupled to a thermal bath, the bare lead-electron Green's func- 
tions are in a voltage driven thermal equilibrium. They have 
the matrix form in the space of electron spin ( [G Q (t, t')] arT i = 

-i((c a a(t) C la' (*')». Where C aa = En C an CT ) 



G„(w) 

Gq1,2M 

h 0a (w) 
p QK (e) 



= / de 



G a i(w) t + G a2 (u) m Q0 ■ <r, (13) 

±[G a+ {u>)±G a -(w)], 
/w(e) 

LU-£ + iO+' 

2iri h 0a (uj) p aK (uj), 
h eq (uj — fj, a ) = — tanhj 
p(e - nB a - fi a ), 



L P - Po 

2T 



where k = ± and p(u>) is the bare electron density of states 
(per spin) which we assume to be similar in both leads. We 
consider here only local Green's functions, so the eigenstate 
index n can be dropped. Note, that for fully polarized elec- 
trons (B a ^> D) we obtain G a i(u>) = Goa(u)), while 
for weakly spin polarized electrons (B a <C D) G a 2{oj) <C 
G al (u). 

The propagators for the Majorana fields can be compactly 
written as G^i, t') = —i^J} l xit)Vv{t')Jl- The bare propaga- 
tors satisfy G xx = G yy , G yx = — G xy , G xz = G° zx = G yz = 
G° zy = 0, due to the symmetry under rotations about the z- 
axis. A useful property holds for the Keldysh components of 
the spin Green's functions: Since Majorana fermion operators 
are Hermitian, from the properties of the commutator, it turns 
out that G I ^ lJ (t 1 1') = —G^Atf, t). In the frequency uj repre- 
sentation, for time-translation invariant solutions, this means 
that G% v (uj) = -G*J—u), and therefore G^(w) is an odd 
function of u>. 



To write down Dyson's equation, it is convenient to use 
a basis in which the bare propagator is diagonal. Defining 
the canonical fermion operators / = {rj x — ir)y)/V2 and 
= (Vx + iT]y)/V2, with the associated Green's functions 
G f p(t) = -i«/(i)/t(0)» and G fff (t) = -*«/t(t)/(0)». 
one has the following expressions for the bare propagators in 
the (/, /' ,i) z ) basis in the steady state in the lowest order of 
perturbation theory 



G°» = 



G%{u) 



is 
1 



G z f (w) 



B 
1 



B 



/t/M 



G 



G 0K 



2ih z (uj)s 
UJ 2 + s 2 ' 
2ihf(to)s 
' (u - B) 2 + t 
2ihf\- (lo)s 
' {oj + B) 2 + i 



(14) 



where s — ► + is a small regulator (a width due to an infinites- 
imally small coupling to a thermal bath), and B is the total 
effective magnetic field acting on the spin (see next section). 
In the thermal equilibrium state h z (u>) = hf(u>) = hp (u) = 
-tanh(w/2T). 

The perturbative calculations of various physical quantities 
are very conveniently performed by using the matrix form of 
the Green's functions in the Keldysh space. Graphical repre- 
sentations can then be done by Feynman diagrams with lines 
representing matrices of Green's functions, where at each ver- 
tex of an internal point is assigned an additional + or — factor 
(due to the opposite direction of time integration for the points 
on the — part of the Schwinger-Keldysh contour). The interac- 
tion vertices of Majorana fermions with electrons follow from 
the form of the tunneling Hamiltonian Q 



rri(ex)c 
1 aj3 a 



1 a/3 3 



T. 



(ex) 



a/3 



T 



(ex) 



^(/ + /-// f ) + -^(/<7+-/t<7-)7fc, 



T. 



(so) 



T 



^/2 

(so) 



^U/ t /-// t ) + ^k/z + -/ t r>7 



V2 



a x ±ia v , I = l x ±iL 



(15) 



accounting for the relation ff< + f'f = 1. We represent the 
electron and spin Green's functions as well as the vertices in 
Fig.0 

The interacting Green's function for the spin, G, satisfies 
the Dyson's equation 



G- 1 (w) = Go 1 (w)-S(w) ) 



(16) 



for the 6x6 matrices (in the tensor product of the x, y, z 
space and the Keldysh space), where the free propagator is 
given by Eqs. dl4> . To lowest order, the imaginary part of 
the self-energy is of order T 2 pPQ (its w-dependent real 
part induces a trivial small energy shift that we neglect in the 
following). Calculating G(u) and neglecting terms of order 
TapPo we can take the 2x2 blocks £//, S/t/t> S/ z , S/t z , 
S z f and E„ft as zeroes. Then, G _1 ((jj) has nonvanishing 



j z f anu z-> z f\ as zeroes, 
off-diagonal blocks ff' and f'f, 



G 



IP 



B 



Ef /t (o;) 



U) 



-Sf /t ( W ) 



(17) 
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FIG. 2: Diagrams for: -the lead electron Green's function (solid line); 
-the spin Green's function; -vertices for the exchange interaction. 



and the diagonal block zz. Inverting such block matrix 
G _1 (cj), we obtain the matrix G(u>) with a similar block 
structure, where the block fp is given as 



-B-T, 



^7? 



-B-T, 



1 



(18) 



while for Gyt/ we replace / — > f< and B — > —B. For G Z2 
we replace zz for /t/ and put B = 0. 

Parcollet and Hooley |y] have found that if one starts from 
the equilibrium distribution functions for Majorana fermions, 
the perturbation theory breaks down when one takes the limit 
of zero coupling to the leads prior to taking the limit of 
zero coupling to the thermal bath. The reason for this non- 
commutativity of the limits is that in a non-equilibrium steady 
state the distribution function of the system can deviate from 
the equilibrium one significantly in the long-time limit even 
for very weak tunneling. Hence, at long times it cannot be 
calculated perturbatively. One can, nevertheless, use the per- 
turbation theory built upon the appropriate bare Green's func- 
tions with correct zeroth-order distribution functions h z (ui), 
and h f * t (w) which are stable with respect to weak perturba- 
tions. These distributions should be obtained self-consistently, 
i.e., h z (u) and hftt(w) computed using second order pertur- 
bation theory have to be the same as the zero order distribution 
function. Hence, assuming that after a long time the system is 
in a steady state, we can define a distribution function hf(u) 
obeying the self-consistency equations 



h f (u) 



h f (u)[Gf r {Lj) -Gf r (u,)], 



yA 



yR 



(19) 



where £<t / (w) is calculated in the second order perturbation 
theory with respect to the electron-spin interaction. hfi(ui) 



FIG. 3: Diagrams for the self-energy due to the electron-spin inter- 
action that contribute to the effective dc magnetic field. 



and h z (ui) have to be obtained in a similar way. Therefore, 
one needs to take the zeroth-order Green's functions for the 
spin in the form of Eq. J14-I . calculate the second order self- 
energies Yip Jul) and S^ t Acj)— £j? t ^(cj), and solve Eq. dl9l l. 
This procedure, as well as taking the limits s — > + and 
T a p —> 0, was carefully outlined by Parcollet and Hooley 

A useful relation between the distribution functions hf(ui) 
and hp (ui) can be derived. Using the definitions in Eq. dl 11 1 
and the commutator and anticommutator properties one can 
show that 



-Gf ft (co), (20) 
Gfi f (-u,)-Gf if (-oj) - Gfp(w) - Gfp(w), 

and using the definition of Eq. ( I19> one obtains that 

h f (uj) = -hp(-w), h z (u) = -h z (-ij), (21) 
so that h z (0) = 0. 



IV. EFFECT OF THE TUNNELING ELECTRONS ON THE 
SPIN DYNAMICS 

A. Effective dc magnetic field acting on the spin 

The diagrams shown in Fig.|5]lead to a renormalization of 
the effective dc magnetic field acting on the spin when the 
electrons are polarized. To first order in the electron-spin in- 
teraction (the first diagram in Fig. [3jl the additional dc mag- 
netic field caused by the tunneling electrons is given by 



(jB t = 2Re 



a,n,cr,<7' 

s m nT( ex ) 



(22) 



neglecting spin-orbit corrections. Here s a is the degree of 
the electron spin polarization in the lead a, i.e., the ratio of 
the electron magnetization to the optimal magnetization. This 
additional magnetic field originates due to the exchange inter- 
action between the spin and the electrons in the leads, and it 
is voltage independent. 

We may estimate this contribution using the expression for 
the tunneling resistance Rt = h/(47re 2 T§pl). For Rt =100 



6 



LA-UR-OO-XXXX 



February 6, 2008 



Submitted to Physical Review B 




(3 f3 
FIG. 4: Diagrams for the self-energy to second order in T a p. 

Mohm we obtain T p w 1.6 ■ 10~ 3 and taking « 0.1T 
we estimate at po ~ eV -1 the additional magnetic field 
Bt ~ 16s Q T. When electrons are polarized by the mag- 
netic (or exchange) field B a such that B a po <C 1, we estimate 
Bt ~ 10 _3 -B a . The second diagram in Fig.[3]is smaller by a 
factor TqPo and, thus, may be neglected. 

In the following we choose the z-axis along the total field 
B = Bo + By- By m Q we denote the direction of elec- 
tron spin polarization in this new coordinate system. In the 
case that rriRo = mLo = mo, and mo has the components 
mo x ,moz in the coordinate system with the z-axis aligned 
along the bare field Bo, we get that m has coordinates m x , m z 
in the coordinate system where the z-axis is aligned along B 
with 



y/B* - (Bo • mo) 



B\ + 2B T B 



(B • m ) 



o 
Bt 



m 



(23) 



B^ 



2_Bt i Bq ■ mo 



Self-energy of Majorana fermions and steady state 
distribution function 



The diagrams shown in Figs.|4]determine the imaginary part 
of the Majorana fermion self-energies which depend upon T, 
V, m Q and B. We start by considering only the interaction 
of the spin with the tunneling electrons, assuming T env = 0. 
In this case, the spin is precessing freely between consecutive 
passings of tunneling electrons which change the phase of pre- 
cession randomly, preserving time-translation invariance on 
the average in the steady state. In the quasi-classical approach 
id these passings are replaced by an effective classical mag- 
netic field representing white noise. On the other hand, in 
the quantum description the electrons are treated as tunneling 
particles with initial energy e„ + V and final energy e n > . The 
tunneling process now has similarities with the scattering of 
an individual particle by the spin. The difference is that both 
the initial and final energy of the tunneling electrons are not 



fixed, in our model being restricted only by the bands of the 
leads-electrons. 

The diagrams shown in Fig.|4]lead to the self-energy 



„_l / n v—v „9 f dudvde „, 



{iTr[a z G+-(u)a z Gp + (v)]G+-(e) + 
2TT[ fJ +Gt-(u) f j-G- p + (v)]Gt~(e)}. (24) 

The expression for ^Jp can be similarly obtained. Note, that 
SjTf depends on hp (— B) via Gp,, while depends on 
hf (B). Then, Eq. dl9l is a closed equation for hf(u>). The 
other diagrams in Fig.|4]give 

/ > ^ f dudvde „. 

a/3 



{Tr[a-G+-(u)a+Gp + (v)]G+-(e) + 
Tr[a+G+-(u)a-Gp + (v)}G+p(e)}. (25) 

The traces for the spin matrices in the expressions for the self- 
energy can be performed with the help of the relations 



7(5 

-m Q • m ( 8^ + m att mp v + m au m 0ll , 
m av Tr[a^,a u ] = 2m all , Tt[a x a y a z ] = 2i.(26) 

V 

The transformation from the "± basis" to the "LO basis" for 
the self-energy is the same as for the inverse Green's function 

which implies that + £±~ + + £7~ = and, 
therefore, S K = -(£+" + £"+) and H A - Y> R = £"+ - 

s+-. 

In this way, for the fully spin-polarized case, m a = mp 
1 . we obtain 



,1*13- fJ, a + u±B 



a,/3 



[(ji a -lif,-uTB)Thf(B) T</>(- 



T 



)] 



-(1 - m az mp z ± m az =F mp x ){np - fx a + lu)} (27) 



and 



a,/3 



(28) 



{(1 + 2m az m (iz - m„ ■ m/j 

, M/3 - Ma + ^ ± S 



±hf(B)(jip -/j, a + ui±B)] 



-(1 - m az mj3 z ± to qz T mp z )T(j)( — )}, 
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FIG. 5: Upper panel displays the steady state distribution function 
hf(B) in the situation m-Rz = m^z — as a function of b and v, 
for two different values of 8. Lower panel shows contour plots of 
the same function. The case v = corresponds to the equilibrium 
distribution h sq (B). 



with the property 



ImE 



A 

//+,/+/ 



-ImS 



R 



-(E A - S fl ) //tj/t/ . 



Assuming that D V,B, T, we calculated the integrals over 
electron energies, 



duG+i {u)G f3 + {u + w), 



in the limit of large electron bandwidth, D — > oo, and us- 
ing the relation 1 ± h eq (uj) = 2np(±Lo) (where np(x) = 
[exp(x) + is the Fermi distribution function) the result- 
ing function cf>(x) is 



4>(x) 



tanh(a;/2) 



(29) 



In the weak polarization limit, B a <C D, there are correc- 
tions to these expressions due to terms consisting of integrals 
of the form 



duG+ 2 {u)G^{u 



w). 



To estimate them we assume that near the Fermi energy the 
electron density of states may be approximated as p(e) w 
po(l + e/D). Then, these integrals are proportional to the 
small parameter (B a /D) 2 and they vanish in the limit D — > 
oo. We obtain the expressions for self-energies in the weak 
polarization limit from the fully polarized results by putting 
m Q = and multiplying the resulting expressions by a factor 
4 that accounts for both spin projections of the electron. 



In the next section we show that the imaginary part of the 
self-energy E f ft (ui), T± = ImE * p (cu), describes dephasing 
of the spin precession caused by the tunneling electrons. It 
determines the current noise power as shown later. 

Therefore, at frequency co = ±B, Eq. dl9l gives us a self- 
consistent equation for hf(B) = —hp(—B) in the case of 
full spin polarization. This equation has a unique solution for 
the distribution function (b = B/T, v = V/T) 



hf(B): 



2b(l - m Rz TO Lz )~2v(m Rz - m Lz ) + b 
<j>+(l - m Rz m Lz )-(j)-(m Rz - m Lz )+0(b) 6' 
^ R (l-mL)+T L 2 L (l 



m iz) 



T 2 

RL 



0(v + b)±0(v-b), 



(30) 



provided that we do not have one of the cases m,Kz = rwLz = 
±1. In the latter two cases, the self-consistency equation is an 
identity and the spin steady state can be any. When m Q = 
we reproduce the results of Ref. |2;] for unpolarized electrons 
(see Fig. [5}. If we account for environment-driven relax- 
ation, r en v 7^ 0, the distribution function will be closer to 
that in thermal equilibrium. When such relaxation dominates, 
the distribution function hf(ui) coincides with the equilibrium 



C. Average spin magnetization and the spin-spin correlation 
function 

We compute first the average spin magnetization (S^(t)) = 
— {i/2)e^ vl {r] v {t)ri 1 {t)) (time-independent in the dc steady 
state), where e^ 1 is the antisymmetric unit tensor. To lowest 
order, the values (S x ) and (S v ) are proportional to the small 
parameters T^p\ and can be neglected. For (S z ), we get (to 
zeroth order) 



(S z 



h f (B). 



(31) 



When rriR, = rxiL (except m az = ±1), we obtain the same re- 
sult as for unpolarized electrons found by Parcollet and Hoo- 
ley 



= ! tax* (H 



#b)(2 + fl) 
4>+ + 4>(b)6~ 



(32) 



Here = (T^ L + Tr R )/T| l . At T = 0, and V < B we 
get (S z ) = 1/2, i.e., the spin is in the ground state because 
electrons do not have enough energy to flip the spin. When 
V > B the electrons can flip the spin, thus reducing the spin 
magnetization to (S z ) = (l + 0/2)/(2x+0), where x = V/B. 
It drops as 1 /V for large V. 

As noticed by Shnirman and Makhlin lUrjll the spin-spin 
correlation function S X x(t) = (S x (t)S x (0)) may be ex- 
pressed via the /-fermion Green's functions. We rewrite the 
spin operators in the following fashion 



S x = -(/ t + /)f/2 = -f(/t + /)/2, 
S y = - f)f/2 = tr(p - f)/2, 
S z = -7?*t/\/2 = -tt) z /V2. 



(33) 
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The operator f = y/2(l-2f^ f)r) z = i2V2r) x r\ y T\ z commutes 
with the spin operators 5 M and hence with the Hamiltonian, 
while f 2 = 1. (Notice that the operator f maps the orthogonal 
subspaces TL\,± «-> Ti.2,± of section ITTT1 among themselves.) 
As a result, 



S xx (t) = (S x (t)S x (0)) 
1 



(34) 



i ([/ t W + /W][/ t (0) + /(0)]) 



■ [G-+ (t) + G-+(t) + Gj+ f1 (t) + G-+(t)}. 



Similarly, 

Syy(t) = 

S zz (t) = 



- i <[/ t (*)-/(*)][/ t (o)-/(o)]> 

\[G- f ^) + G-+{t)-G-+ p {t)- 



(35) 



and, in general, S^ u (t) = \G^(t). 

To zeroth order in perturbation theory we obtain 

S xx {uj) = S yy (uj) = - \G-±(u) + G7 t +H] 



(36) 



= -{[1 - h f (B)]5(u> -B) + [l + h f (B)]5(oj + B)}, 

while S zz {uj) = (tt/2) 8{uj). 

To describe the contribution which displays a peak at the 
Larmor frequency in the correlation function ^^(u;) we need 
to account for the imaginary part of the self-energy. For 



G 



//+,/*/ 



we have 



lu - B - £ 



Ft , A ' 



(37) 



B - S 



R,A 



At the weak tunneling condition, |Trl| 2 Po ^ 1' wnen the 
relaxation due to the environment is negligible, we can also 
neglect the contributions from Gff and Gyt/t to S xx with 
respect to those from G fp and Gp f. We obtain 



(l-h f (B))T ± , (l + h f (B))T ± 



{u-Bf+Y\ {u + Bf 



(38) 



where, at T = 0, for fully polarized electrons along the x-axis 
we get 

T X (B,V) = lmEj if (B) 

= jT* L p 2 [\V + B\ + \V-B\+9B}. (39) 

When V < B this gives the width of the precession resonance 
due to quantum fluctuations of the current (the current oper- 
ator does not commute with the Hamiltonian, and the current 
fluctuates between the leads and the spin), 



Tf{B) = \T^plB{2 + e), 



(40) 




FIG. 6: A diagram for the current-current correlation function that 
does not contribute to the peak at the renormalized Larmor frequency 
but gives a spin-dependent correction to the shot noise. 



while for V > B the decoherence rate has an additional con- 
tribution due to the voltage induced transport current 

T ± (B, V) = T [ l\B) + -T^fiiV - B)Q(V - B). (41) 

We estimate Tj_ « 10~ 5 F, in a junction with T RL = 0.1T 
and tunneling resistance Rt = 100 Mohm. To determine T± 
in the unpolarized case we need to put m Q = and multiply 
by 4 in Eq. (39). 

Thus, we see that the spin-spin correlation functions 
S xx (u>) and S yy (uj) exhibit oscillations at a renormalized Lar- 
mor frequency ujl, with a peak width depending on V, B and 
T. When the relaxation due to the environment dominates, 
B 3> r env ^> T±, we need to replace T± by r env in Eq. ( l38l . 
The important point is that T± provides a lower bound for the 
decoherence rate T, i.e., for the width of the peak at lo^ in the 
spin-spin correlation function. 



V. EFFECT OF THE SPIN ON THE TUNNELING 
CURRENT 

A. The average current 

Let us now turn to the evaluation of the average current and 
current power spectrum. The current operator is given as 

I(t) = -ieJ2 [cL CT ( i )(^L)-' c Wa'(*)-H.c.]. (42) 

n.n' ,{7 : a' 

The average current (involving the complete spin-electrons 
Hamiltonian H) is determined from 

I = (I(t)) = 2e Im (^W^rO-'Cl^)). 

n,n' ,<r,cr' 

(43) 

In the following, we will assume that Ht is turned on adiabat- 
ically and we are interested only in the average current up to 
second order in TCt- (Remember that we are only interested 
in the dc steady state solution.) The following derivation is 
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similar to that explained in Ref. 11511 for the quantum dot tun- 
neling. Denoting 



Mt)= J2 c L ct (*)(^rlW'c l „ v ,(*), 



(44) 



we obtain (assuming that the leads are not superconducting) 

(I(V)) = -2elm{U B (-V)}, (45) 
U R {uo) = -i dte lujt 0(t) <[i(i),it( )]). 



Assuming the dependence p(e) = po near the Fermi energy 
we obtain for fully polarized electrons, B a D, 



I(V) = I Q (V)+I S (V)-{S), 
I (V) = 7re(l + m R • m L ) T 2 p 2 V, 

I?>{V) = 27re(m R + m h )T T R ^p 2 V. 



(46) 



For weakly polarized electrons, B a = B -C D, approxi- 
mating p(e) « po(l + e/D) near the Fermi energy, we get 

|I S | cx eT T^4(B/D)V, while I = 4TreT$p>V. 

.(ex) 



When (m R , hil) J- B the term of order TqT^£' vanishes 
because (S) is parallel to B. The next spin-dependent term in 

the average current, 1^, is of order (Tj^) 2 and, in the fully 
polarized case along the 2>axis, is given by 



e) = 



2e(T^) 2 Re{ 



{G R (t' - t)G L (t - t') S xx (t - t') - 
G R (t' -t)Gt(t-t')S xx (t' -t)}}. 



(47) 



Here G a (t) is the causal Green's function. Using Fourier 
transforms and Eq. (I36> the result, at T = 0, is 



(2) 



hf(B)B] Q(V — B). (48) 



Hence, when V < B the average current is not affected by 
the spin because electrons have not enough energy to flip it. 
The spin remains in the ground state and thus cannot affect 
tunneling electrons and be probed. In the case that there is 
no reflection, i.e., 9 = 0, the spin-dependent contribution is 
related to the spin-flip rate IV 



(49) 



because the spin-dependent contribution to the current is de- 
termined by the matrix element of the operator S x between 
states with opposite spin projections. Hence, each tunneling 
electron needs to flip the localized spin in order to go through. 
In the large V limit we obtain T|| — ► T± in agreement with 
the quasiclassical results llj|3j]. 

(2) 

At T = 0, the contribution leads to a discontinuous 
derivative of I(V) at V = B due to the opening of a new chan- 
nel corresponding to the tunneling of electrons via the spin. 



The change is more pronounced in the second order deriva- 
tive of I(V) with respect to V. Accounting for broadening of 
the spin energy, Tj_, we derive 



d 2 I 
dV 2 



2e( r RL X) )Vo E S^{V){m R4i + m^) 2 . (50) 



H=x,y 



We see that this second order derivative of the average tun- 
neling current shows a peak (see Eq. i3bt ) at V = B. Hence, 
measurements of the I-V tunneling characteristics provide in- 
formation on the spectral density of the spin-spin correlation 
function. Similar behavior for the average current with re- 
spect to bias voltage was used previously for inelastic tunnel- 
ing spectroscopy of phonons 11711 . It is based on an equation 



similar to Eq. (I50i relating d 2 I/dV 2 to the phonon density of 
states. 

Note, that in the case of unpolarized electrons I s is small, 
but the second derivative of I(V) also shows a similar peak 
fjlil . However, now its amplitude has an additional small fac- 
tor T 2 p 2 . 



B. The shot noise power 

Next we calculate the current-current correlation function. 
The particular fluctuation operator one needs to evaluate de- 
pends upon the nature of the measurement itself fil^l . see also 
12011 . In the general case the result of measurements depends 
on the correlation functions describing emission, (+), and ab- 
sorption, (-), by the tunneling contact 



We obtain 



dt e 



V{t), V{t) = (I(0)I(t)}. (51) 



P±(u)=Ps±(u)+V p± (u), (52) 
V s (t)=e 2 £ [((4,„ ff (0)c R „, 7 (t)) e (f RL )^(0) 

n,n / ,m,m / ,(T,cr / ,7,7' 

x(c Lmff ,(0)c[ m , 7 ,(t)) e (f RL ) y7 (i)) s + (R<-> L)], 

where the term V s (t) describes the shot noise to second order 
in the matrix elements To and Trl, while V p (t) represents the 
contribution which displays a peak at the renormalized Lar- 
mor frequency of order (ToT R l) 2 (see next subsection). We 
calculate Vs (t) in the framework of standard time-dependent 
perturbation theory (as we did to calculate the average cur- 
rent). In the equation above, (...) e means average over elec- 
tron degrees of freedom, after the perturbation due to the spin- 
electron interaction is accounted for, while (...)« means aver- 
age over the localized spin degrees of freedom. Averages over 
spin operators are performed using the known spin correlation 
functions determined in previous sections. 

Contributions proportional to Tq and To Trl gi ye the stan- 
dard expression 12111 modified by the presence of the localized 
spin, in the same way as the average current Eq. d46t . For fully 
polarized electrons along the .x-axis we obtain, at T = 0, 



Vi°)(u) = 2ire 2 (T 2 + 2T T nL )p 2 \V - u\ 



(53) 
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For unpolarized electrons we need to multiply this result by 
a factor 4. We obtain for the symmetrized shot noise power 

(lo < V) 

[P^+V^]/2 = e(I +lP). (54) 

The terms proportional to T^ L describe the effect of the spin 
dynamics on the shot noise, Vs . A diagram contributing to 
Vs is shown in Fig. |6] For electrons fully polarized along 
the x-direction we derive for the emission part at T = and 

r± = 

7>&>H = \e\T^fplQ{V-B)[{l-h f {B)) (55) 
xF(-w + V - B) + (1 + h f (B))F(-w + V + B)}, 

where 

F(x)=x@(x). (56) 

Emission is possible at frequencies uj < V — B, when the 
spin is in the ground state, and frequencies lu < V + B, when 
the spin is in an excited state. The second order derivative 
of Vgl (u>) with respect to ui, when V > B, has peaks with 

amplitudes 7re 2 (T^ c L x) ) 2 po(l ± h f {B))/2 and widths T x at 
frequencies u> = V ± B. This provides a way to detect the 
presence of the localized spin. 

The absorption part at T = and T± = is given as 

0{V - B) [(1 - h f {B)){uj + V - B + F(lu -V - B)) 
+ (1 + h f (B))((j + V + B + F(uj-V + B))} 
+29(5 - V) [F(lu + V - B) + F{uj -V - B)] }, (57) 

with the result 

[V [ £+V^]/2 = elfl (58) 

for uj < \V — B\. Figure^displays the emission and absorp- 
tion contributions. 

For unpolarized electrons we derive 

viT' 1] H = #>S(«) + ^{T^fpl (V - B), (59) 

where V^(lu) is given by Eq. d55l > with the same hf(B) as 
for fully polarized electrons. This noise power also exhibits 

anomalies at frequencies u> = V ± B. We see that Va (u>) 
has no peak at the Larmor frequency neither for polarized nor 
for unpolarized electrons contrary to the results of Ref. Jl4ll . 

The fourth order terms, Tq 7r L , m T^i 10 ) correspond to 10 
skeleton diagrams some of which are shown in Figs. [8] and 
|9] In these 4- vertex diagrams the external frequency (cu) (via 
the current vertex) may enter in any of the two vertices. Di- 
agrams in Fig. [S](all six), as well as the diagram shown in 
Fig-El contain integration over the energy e entering in all the 
electron and spin lines of the diagram. Hence, these contribu- 
tions change at least on the frequency scale of order V, B as 
the term Vs^ioj)- They are small in comparison to Vs^iuj), 
and may be neglected. 




1 2 3 4 5 

UJ [units of B] 




X=0 12 3 4 5 

UJ [units of B] 



FIG. 7: Emission and absorption contributions to the shot noise 
power to order Ijf L , at T = 0, for similar electron polarizations in 

the leads m 1 B. T>i± is in units (7rBe 2 (T| L po/4). For x < 1 the 
emission contribution vanishes. Notice the different regimes, charac- 
terized by a change in the slope, that appear because of the presence 
of the localized spin. 



(•) and (■) vertices for current 




FIG. 8: Two of the six diagrams for the current-current correlation 
function that do not contribute to the peak at the renormalized Lar- 
mor frequency to order Tq Irl- The external vertices can be any two 
of the four vertices in the diagram. 



C. The peak at the renormalized Larmor frequency 

To derive this peak we consider the 4 diagrams of the type 
shown in Fig. [9] which give contributions of order TqT^ l . In 
these diagrams integration over energy e in the spin lines is 
independent of the integrations in the electron lines and they 
provide a contribution to the current-current correlation func- 
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tion V p (u>) which has the same peaks at u> = ±B as the spin- 
spin correlation function. To calculate this contribution we use 
the Keldysh technique as done by Shnirman et al. fl^l . We 
will assume that to study the peak at the renormalized Larmor 
frequency experimentally one can measure the instantaneous 
values of the current over a long period of time and then obtain 
the symmetrized current power spectrum defined as 



V{lj) = -[V + (u)+V-{uj)]. 



(60) 



To calculate the first term in Eq. \5\\ . we put t = t- and 
= 0+. For electrons fully polarized along the x-axis the 
first term is given as 

(J(t)/(0)) = (-T 2 S(t)S(0) -ToT^A^S^BiO) 
+B(t)A(0)S x (0)} - T^A(t)S x (t)B(0)S x (0)), (61) 

£(*) = £4n*(*)<W*)> 

n 

A(t) =^2c+ nx (t)a x c Lnx (t). 

Now, in the framework of the Keldysh technique, we expand 
B(t)B(0) up to (T RL ) 2 order in the first term, up to T Tr. l 
in the second and third terms and up to Tq in the last term. 
Let us consider the first term which corresponds to the second 
diagram in Fig.|9] This diagram leads to the contribution 

-(T T RL ) 2 E A r £(- w ) A ua^>)A'+H. 

/X,V=zt 

Here the product pv is the sign factor (+1) if fx = v and (—1) 
otherwise, while 

A RL H = J ^[G{r(«( e - W )], (62) 
A^>) = [A R >)-A£ R (o;)], (63) 

where the trace is taken over spin variables. Further, we de- 
note by (u>) the combination of Majorana fermion Green's 
functions which are represented by broken lines in the dia- 
grams in Fig. [9] 



= J ^G%(e + u)G%(e) = (64) 
J ^G%{e + u;)G%(e) = J ^G^G^e-u;), 
with the property 

S'"'(u)=S vl *(-u). (65) 

Adding the term with interchanged R and L and symmetrizing 
with respect to time we obtain for the first term 

P P W M = -i(ToT RL ) 2 E^[A^(^)5^(c)A^(c) 



Similarly we calculate the other three terms 

1 

2 



(T T RL ) 2 E HA^(0)<S-»A"+H 



/ 1 ■ V 



+A- f *(0)S- u (-uj)A» + (-uj)}, (67) 
Vf\u>) = -i(T T RL ) 2 E^[A^^)^-(c)A^(0) 



/ 1 . V 

+A-^(-uj)S^-(-uj)A"+(0)1 

^ 4) H = -i(T T RL ) 2 E^A^(0)[5" + H + 

fJ,,V 

S-+(-uj)}A v+ (0). 



(68) 



(69) 



Let us calculate S^foj) using the expressions for G + and 
G H in terms of the zeroth-order functions G R , G A and G K 



1 f de 

4 J ^ C vi 



{G&(e + w)G£(e)+ (70) 



(e + w)G£(e) + G* (e + w)Gft (e) + 
Gf z (e + W )[G*(e)-G£,(e)] + 
G*(e + u,)[G^(e)-G*(e)]}. 

Since /i 2 (0) = 0, we have G^ z (ui) = 0. Therefore, the contri- 
bution to S '"(cj) from the terms containing G^j is negligible, 
but the last term does contribute. We get 

5-+H =5 +_ (-w) (71) 

= £{[!- ft/MW" - B) + [1 + + B)}, 



5~(w) = -J%(w)P 



B 2 



(72) 



(73) 



where S + (uj) = 6(u> — B) + S(uj + B). Next, we take into 
account the broadening of the resonances at w = ±B by re- 
placing 



tt8{uj ± B) 



(w ± Bf + r 2 ± 



(74) 



+A-^(-w)^ i/ (-w)A I/ +(-w)], 



(66) 



Such procedure follows from the fact that accounting for 
higher order terms of perturbation theory in T R l we need 
to replace the product G%" (t)G%% (-<) in Eq. (|64j by 
-(Vy{tfi,)Vy{0v)Vz(0v)r]z(tfj,)}. This function has the same 
frequency dependence as Gfp (t) because from one hand 

(S x (t)S x (Q)) = -(Vy(t)v*(t)riy(0)ri z (0)}, (75) 

while from the other this spin-spin correlation function is de- 
termined by Eq. ( 134k At T = we obtain, neglecting terms 
of order uj / D, V/D, 

A-+(w, V) = A-+(u,, — V) = 27r(c + V)Q(u> + V)p 2 , 

-2tt((j + V)Q(-u - V)pl, 

A++(u;,V) = A++(lj,-V) = A^(u,V) = 

AI k {oj,-V)=tx\oj + V\pI (76) 
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FIG. 9: Two of the four diagrams contributing to the peak at the 
renormalized Larmor frequency in the current-current correlation 
function to order TqT^ l . There are two other diagrams with dif- 
ferent current vertices. 



Summing up all terms for P p (u>) we obtain zero for the am- 
plitude of the peak at V < B. This same result was obtained 
by Shnirman et al. l22ll . It means that quantum fluctuations of 
the spin in the ground state which are present in the spin-spin 
correlation function cannot be probed by tunneling electrons. 

When V > B we obtain, taking into account the broaden- 
ing of the resonance, 

V p (cj, V, B) = P(V, B) V± (77) 
1 j_ + [uj — a)- 

P(V, B) = 4(irT T RL p 2 ) 2 [V 2 + VBh f (B)}. 

More generally, the coefficient 4 is substituted by [m/jj_ + 
m^] 2 for fully polarized electrons in the leads, while for 
weakly polarized electrons we get the factor |mp{^SR + 
m L± s L| 2 - Accounting also for spin-orbit interaction we ob- 
tain an additional term in P(V, B) which is proportional to 
the small factor [{2B a /D] 2 + /3 A (t 2 + t 2 )}. 

We see that P p {w) exhibits a peak at the renormalized Lar- 
mor frequency with a width T± only when the electrons are 
polarized in a direction which is not parallel to Bo- This peak 
is caused by precession of the spin excited by tunneling elec- 
trons. Tunneling of electrons excites the spin, resulting in the 
modulation of the tunneling current, but they also cause deco- 
herence of the spin precession with a rate T±. 

Therefore, the spin precession cannot be seen in the current- 
current correlation function without spin-orbit interaction 
when the leads-electron polarization is parallel to Bo, or the 
electrons are unpolarized, despite the fact that the spin-spin 
correlation functions for the S x and S y spin components ex- 
hibit oscillations at u>l. The reason for this is that to observe 
oscillations in P(uj) one needs to probe with the current a 
component of the spin perpendicular to B. For that the elec- 
trons need to have a component of polarization rriR or rriL 
that couples to S x or S y . Spin-orbit interaction leads to a sig- 
nal when the electrons are unpolarized or Bo |j m a o if Bo 
is not parallel to 1. However, the signal in that case is much 
weaker than in the case of fully polarized electrons. 



30 




2 4 6 8 10 



x = V/B 

FIG. 10: Width of the peak as a function of voltage for the case 
where electrons are fully polarized in the direction perpendicular to 
the effective magnetic field acting on the spin. 

D. The signal-to-noise ratio R 

For the symmetrized signal the signal-to-noise ratio reaches 
its maximum at low temperatures when r env <C T±, while 
m.n || rri£, and both are perpendicular to B. Indeed, elec- 
trons with components of the spin polarization parallel to B 
do not contribute to the signal but they do enhance T±. When 
mR± = mh± = 1 the relaxation rate T±(V 7 B) (neglecting 
relaxation due to the environment) is determined by Eq. i39i . 
The signal-to-noise ratio at the peak position, for the case of 
perpendicular spin polarization at T = and V > B, is given 
by the expression 

= P P (B) = P(V,B) = V + h f (B)B 
2nVT 2 p 2 r ± (V,B) V + 9B/2 ' 

(78) 

Here P p (B) is the height of the peak at the renormalized Lar- 
mor frequency in the symmetrized current-current correlation 
function, while pf^ is the current power spectrum for the shot 
noise neglecting the smaller contribution Vs . 

R reaches its maximum value i? max = 4 for V ^ B. In the 
case of symmetrical leads, Trr = Tll = ?rl we get 9 = 2, 
and the function R(V/B) for = 0, 2 is shown in Fig. II II 

When the spin relaxation due to the environment becomes 
dominant, r env ^> T±, the signal-to-noise ratio for the case of 
fully polarized electrons is smaller by the factor r^/r en v <C 
1. If the perpendicular components of m Q were absent, we see 
that R becomes as small as /3 4 , while for weak polarization in 
the direction perpendicular to B it is as small as (B a /D) 2 . 

VI. DISCUSSION 

We start summarizing our main results for the model of a lo- 
calized spin interacting with the tunneling electrons only via 
the exchange interaction, i.e., when the tunneling matrix el- 
ement contains the term T^cr ■ S in addition to the spin- 
independent term Tq. 
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x = V/B 

FIG. 11: Signal-to-noise ratio -R(x) for fully polarized electrons. 




1 2 3 4 5 6 7 



x = V/B 



FIG. 12: Comparison between the dependencies of T±, I — Iq and 
R vs. x = V/B at T = and 8 — 0, explaining why R reaches 
approaches a maximum at large V. 



1. The Larmor frequency is renormalized from the one 
determined by the external magnetic field Bo to 5 = 
yj Bq + (see Eq. (I22H . when tunneling electrons are po- 
larized. An estimate for Bt is given after Eq. \22\ . 

2. In the steady state the I-V characteristics, namely 
d I(V)/dV , exhibits a peak at V = B when the electrons 
are polarized in a direction different from B. This peak is due 
to the opening of a new channel for electron tunneling via the 
localized spin (spin-flip tunneling). The width of this peak is 
determined by the spin decoherence rate T±, see Eq. J39l>. 

3. When the electrons are polarized and V > B, on the 
background of frequency-independent shot noise caused by 
the discrete nature of the electrons, there is a peak in the cur- 
rent noise power [V+(u))+V-(oj)]/2 at the renormalized Lar- 
mor frequency, caused by the precession of the spin excited 
by the tunneling electrons. This peak is absent when V < B 
because the spin is in the ground state since electrons do not 
have enough energy to flip it. 

4. The signal-to-noise ratio for this peak, R = 
V p (ujl) /vf^ , strongly depends upon the degree of spin polar- 
ization of the leads-electrons and the orientation of this polar- 



ization with respect to the external magnetic field Bo acting 
on the localized spin. The maximum value of R is reached 
when the electrons are polarized in a direction perpendicular 
to the effective magnetic field acting on the spin because in 
this case they probe in an optimal way the precessing spin 
components. When the relaxation of the localized spin with 
the environment is weak, R reaches the values of order unity 
only if the leads are almost fully polarized in a direction close 
to perpendicular to B. 

5. In the case of full polarization and hiro = rriLo -L 
B at low temperatures T <C 5 e ff , V, the dependence of R 
on x = V/B is given by Eq. J78L where 8 comes from the 
leads-electron reflection via the localized spin (see Eqs. ( I30H . 
The maximum value, i? max = 4, is reached when V ^> B, 
see Figs. ^2 an dEl The signal-to-noise ratio is reduced by 
the presence of the electron reflection couplings (6 ^ 0) as 
expected, since this process does not contribute to the current 
but does affect the width of the peak. 

6. When the spin relaxation due to the environment is neg- 
ligible, the peak width in the noise power spectrum (at the 
renormalized Larmor frequency) depends weakly on the ori- 
entation of the polarization. It increases linearly with V, when 
V B 3> T and the tunneling electrons cause spin flips, see 
Fig. [10| At high temperatures T > V, B, the width of the 
peak is determined by T (T± oc T). 

7. In points 2 through 4 we assumed that the spin relaxation 
with the environment, r env , was negligible in comparison to 
the relaxation due to the leads-electrons, r_i_. This assumption 
is crucial to obtain a maximum R of order unity because in 
this case T± oc \T^\ 2 pl, and V V {B) oc {T^fpl/V^, so 

that R does not depend on the small parameter (T^) 2 pQ, 
see Eq. J78I . If the spin relaxation rate due to the environment 
dominates we get R oc T±/T env <C 1. 

8. For weak polarization of the electrons, R becomes pro- 
portional to the degree of polarization, R ~ (B a /D) 2 , where 
B a is the magnetic or exchange field acting on the leads- 
electrons and D represents their bandwidth. This is because 
only polarized tunneling electrons contribute to the signal at 
the Larmor frequency, while all tunneling electrons contribute 
to the broadening of this signal. We estimate B a /D ~ 10 -5 
for B a = 100 G and D = 0.1 eV and this estimate leads to 
R ~ (B a /D) 2 < 10- 8 . 

9. Spin-orbit (i.e., relativistic 1 ) corrections of any type 
IH can only provide a very small signal-to-noise ratio, 
R ~ /3 4 ~ 10~ 8 , in the situation where the exchange cou- 
pling leads to a vanishing signal (i.e., the electrons are un- 
polarized or they are polarized parallel to Bo). The signal in 
this case is determined by the relativistic coupling between the 
spin and the tunneling electrons, while the broadening is de- 
termined by the dominant non-relativistic exchange coupling. 
Hence, for a tunneling process via a single localized spin in 
the steady state, the signal-to-noise ratio becomes of order 
unity only when 

a) the spin is decoupled from the environment, and 

b) the electrons are almost fully polarized in a direction per- 
pendicular to the effective field acting on the localized spin. 

10. We derived the correction to the standard shot noise 
power caused by the presence of the localized spin when the 
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leads-electrons are unpolarized. This correction is small when 
Irl <C To, but it changes behavior at V — B, see Eqs. \55\ . 
and J59I . On the other hand, the standard shot noise power 
increases linearly with V. This new result can be observed, in 
principle, experimentally. 

We are now in position to compare our results to the ex- 
perimental data obtained by Manassen et al. ]5j] and Durkan 
and Welland |@], and decide whether the model of tunneling 
via a single localized spin can explain these data. In these ex- 
periments the polarization of the electrodes was weak, R was 
of order unity or larger almost independently of the orienta- 
tion of the applied magnetic field, and the position of the peak 
in the noise power spectrum V{u>) was found at the Larmor 
frequency corresponding to the applied magnetic field. Our 
results show that a combination of weak polarization and a 
signal-to-noise ratio of order unity are incompatible for the dc 
steady state described by our model, which accounts for both 
exchange and spin-orbit coupling of the tunneling electrons 
and the localized spin. 

Next, at high voltages we reproduced the quasi-classical re- 
sults GHH which are valid in the limit V 3> B, 

From the point of view of the theory of quantum measure- 
ment we see that after a long time, any measurement usually 
leads to a steady state where information on the initial state 
of the quantum system is probably lost. Hence, an important 
question is what is the transient time for a given measurement. 

We note that d 2 I/dV 2 and the contribution V p (u>) in the 
current-current correlation function are proportional to the 
spin-spin correlation function S XX (V) and S xx (cu), respec- 
tively (if the electrons are spin polarized along the x-axis). 
Suppose that one wants to probe a single spin, say A, which is 
part of an ensemble of other spins interacting with each other. 
Measurement of the I-V characteristics, namely d 2 I/dV 2 , 
and current noise power, V(u>), at low voltages will carry in- 
formation on the dynamics of that particular (coupled) spin A, 
whenever tunneling occurs via that single spin. Effectively, 



because of the coupling, that measurement provides informa- 
tion on the dynamics of the whole system and Eqs. (I50i and 
(I77> form the background for such tunneling spectroscopy of 
quantum systems. Note, however, that the subject of study, the 
spin-spin correlation function, is affected by tunneling mea- 
surements (see the way the bare spin-spin correlation func- 
tion, Eq. J36I . was modified in the case of a single-spin sys- 
tem, Eq. d38» . We have shown how this modification may be 
accounted for in the case of steady-state measurements of a 
single-spin system. In a similar way it may be accounted for 
in more complicated systems, and corrections to obtain the 
bare correlation function of the system studied can, in princi- 
ple, be made. 

We note that we presented here a theoretical description of 
the tunneling spectroscopy of a single localized spin- 1/2 sys- 
tem. However, such a description is also valid for any two- 
level system (quantum dot), and Eq. (0 describes the mapping 
between the two-level system and the spin- 1/2 system. 

In conclusion, we described the characteristics of the tun- 
neling current via an isolated single spin 1/2 (two-level sys- 
tem) in the steady state, i.e., in the long time limit after switch- 
ing on the voltage or the tunneling matrix elements. We found 
optimal conditions when the tunneling current carries maxi- 
mum information about dynamic quantum fluctuations of the 
localized spin. We showed what type of information about 
the isolated spin may be extracted from measurements of the 
tunneling current in the steady state. 
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